demo daphnia instructions

# chunk #1
library(ReacTran)
library(simecol)
source("../R/daphnia_deb_phyto.R")
source("../R/daphnia_deb_equations.R")
source("../R/daphnia_deb_lifeequations.R")
source("../R/daphnia_deb.R")
source("../R/daphnia_deb_ibm.R")
source("../R/deb_ibm.R")

time= 1 ref. Vol= 10 abu= 10 

time= 2 ref. Vol= 20 abu= 9.2 

time= 3 ref. Vol= 20 abu= 8.85 

time= 4 ref. Vol= 20 abu= 8.1 

time= 5 ref. Vol= 20 abu= 7.8 

time= 6 ref. Vol= 2.083333 abu= 49.92 

time= 7 ref. Vol= 4.166667 abu= 47.04 

time= 8 ref. Vol= 4.166667 abu= 44.64 

time= 9 ref. Vol= 4.166667 abu= 42 

time= 10 ref. Vol= 0.6127451 abu= 168.096 

time= 11 ref. Vol= 1.22549 abu= 153.408 

time= 12 ref. Vol= 1.22549 abu= 144.432 

time= 13 ref. Vol= 1.22549 abu= 136.272 

time= 14 ref. Vol= 0.4589851 abu= 235.3018 

time= 15 ref. Vol= 0.9179702 abu= 215.6933 

time= 16 ref. Vol= 0.9179702 abu= 162.3146 

time= 17 ref. Vol= 0.9179702 abu= 124.187 

time= 18 ref. Vol= 0.9179702 abu= 108.936 

time= 19 ref. Vol= 2.753911 abu= 70.8084 

time= 20 ref. Vol= 2.753911 abu= 39.21696 

time= 21 ref. Vol= 5.507821 abu= 26.14464 

time= 22 ref. Vol= 11.01564 abu= 14.16168 

time= 23 ref. Vol= 22.03128 abu= 8.98722 

time= 24 ref. Vol= 22.03128 abu= 5.08368 

time= 25 ref. Vol= 66.09385 abu= 2.85957 

time= 26 ref. Vol= 66.09385 abu= 1.57352 

time= 27 ref. Vol= 198.2816 abu= 0.93806 

time= 28 ref. Vol= 198.2816 abu= 0.57494 

time= 29 ref. Vol= 396.5631 abu= 0.39338 

time= 30 ref. Vol= 793.1262 abu= 0.2219067 

time= 31 ref. Vol= 793.1262 abu= 0.1386917 

time= 32 ref. Vol= 1586.252 abu= 0.08573667 

time= 33 ref. Vol= 3172.505 abu= 0.05043333 

time= 34 ref. Vol= 3172.505 abu= 0.04854208 

time= 35 ref. Vol= 3172.505 abu= 0.04665083 

time= 36 ref. Vol= 3172.505 abu= 0.04539 

time= 37 ref. Vol= 3172.505 abu= 0.04318354 

time= 38 ref. Vol= 3172.505 abu= 0.04129229 

time= 39 ref. Vol= 3172.505 abu= 0.03877062 

time= 40 ref. Vol= 3172.505 abu= 0.03687937 

time= 41 ref. Vol= 3172.505 abu= 0.03435771 

time= 42 ref. Vol= 3172.505 abu= 0.03246646 

time= 43 ref. Vol= 6345.01 abu= 0.03026 

time= 44 ref. Vol= 6345.01 abu= 0.02821115 

time= 45 ref. Vol= 6345.01 abu= 0.02679271 

time= 46 ref. Vol= 6345.01 abu= 0.02553187 

time= 47 ref. Vol= 369.7558 abu= 0.2731532 

time= 48 ref. Vol= 739.5116 abu= 0.2569263 

time= 49 ref. Vol= 739.5116 abu= 0.2474606 

time= 50 ref. Vol= 739.5116 abu= 0.2366427 

time= 51 ref. Vol= 331.6196 abu= 0.2834573 

time= 52 ref. Vol= 663.2391 abu= 0.2593333 

time= 53 ref. Vol= 663.2391 abu= 0.2457635 

time= 54 ref. Vol= 663.2391 abu= 0.2352093 

time= 55 ref. Vol= 292.1758 abu= 0.3627952 

time= 56 ref. Vol= 292.1758 abu= 0.3525274 

time= 57 ref. Vol= 584.3517 abu= 0.3080337 

time= 58 ref. Vol= 584.3517 abu= 0.3028998 

time= 59 ref. Vol= 584.3517 abu= 0.2874981 

time= 60 ref. Vol= 62.23127 abu= 1.783669 

time= 61 ref. Vol= 62.23127 abu= 1.671185 

time= 62 ref. Vol= 62.23127 abu= 1.606909 

time= 63 ref. Vol= 124.4625 abu= 1.494425 

time= 64 ref. Vol= 27.47518 abu= 3.748839 

time= 65 ref. Vol= 54.95035 abu= 3.566856 

time= 66 ref. Vol= 54.95035 abu= 3.403072 

time= 67 ref. Vol= 54.95035 abu= 3.239288 

time= 68 ref. Vol= 19.55529 abu= 5.420529 

time= 69 ref. Vol= 19.55529 abu= 5.164844 

time= 70 ref. Vol= 39.11057 abu= 4.909159 

time= 71 ref. Vol= 39.11057 abu= 4.602336 

time= 72 ref. Vol= 15.21812 abu= 6.636825 

time= 73 ref. Vol= 6.503471 abu= 14.91511 

time= 74 ref. Vol= 13.00694 abu= 14.60759 

time= 75 ref. Vol= 13.00694 abu= 14.14629 

time= 76 ref. Vol= 5.993982 abu= 16.1829 

time= 77 ref. Vol= 2.814076 abu= 31.98208 

time= 78 ref. Vol= 5.628152 abu= 29.84994 

time= 79 ref. Vol= 5.628152 abu= 28.78387 

time= 80 ref. Vol= 5.628152 abu= 27.36244 

time= 81 ref. Vol= 1.803895 abu= 65.41401 

time= 82 ref. Vol= 1.803895 abu= 63.19659 

time= 83 ref. Vol= 1.803895 abu= 59.31609 

time= 84 ref. Vol= 1.803895 abu= 58.20738 

time= 85 ref. Vol= 1.803895 abu= 95.34923 

time= 86 ref. Vol= 0.5352804 abu= 177.4771 

time= 87 ref. Vol= 1.070561 abu= 162.5316 

time= 88 ref. Vol= 1.070561 abu= 154.1248 

time= 89 ref. Vol= 1.070561 abu= 144.7839 

time= 90 ref. Vol= 1.070561 abu= 132.6408 

time= 91 ref. Vol= 2.141122 abu= 88.73854 

time= 92 ref. Vol= 2.141122 abu= 64.91925 

time= 93 ref. Vol= 4.282243 abu= 41.567 

time= 94 ref. Vol= 4.282243 abu= 29.89088 

time= 95 ref. Vol= 8.564486 abu= 19.14884 

time= 96 ref. Vol= 8.564486 abu= 12.25993 

time= 97 ref. Vol= 25.69346 abu= 7.472719 

time= 98 ref. Vol= 25.69346 abu= 4.787211 

time= 99 ref. Vol= 51.38692 abu= 2.88011 

time= 100 ref. Vol= 102.7738 abu= 1.517896 

time= 101 ref. Vol= 205.5477 abu= 0.9243598 

time= 102 ref. Vol= 205.5477 abu= 0.5302906 

time= 103 ref. Vol= 616.643 abu= 0.2773079 

time= 104 ref. Vol= 616.643 abu= 0.1670334 

time= 105 ref. Vol= 1233.286 abu= 0.15406 

time= 106 ref. Vol= 1233.286 abu= 0.143519 

time= 107 ref. Vol= 1233.286 abu= 0.1362214 

time= 108 ref. Vol= 1233.286 abu= 0.1329781 

time= 109 ref. Vol= 1233.286 abu= 0.1248697 

time= 110 ref. Vol= 1233.286 abu= 0.1183829 

time= 111 ref. Vol= 1233.286 abu= 0.1151396 

time= 112 ref. Vol= 1233.286 abu= 0.1086528 

time= 113 ref. Vol= 1233.286 abu= 0.09811187 

time= 114 ref. Vol= 1233.286 abu= 0.09081429 

time= 115 ref. Vol= 1233.286 abu= 0.08594924 

time= 116 ref. Vol= 2466.572 abu= 0.08027335 

time= 117 ref. Vol= 1209.104 abu= 0.09097646 

time= 118 ref. Vol= 145.325 abu= 0.6261827 

time= 119 ref. Vol= 290.65 abu= 0.6124205 

time= 120 ref. Vol= 290.65 abu= 0.5780148 

time= 121 ref. Vol= 290.65 abu= 0.5573714 

time= 122 ref. Vol= 74.90979 abu= 1.428385 

time= 123 ref. Vol= 74.90979 abu= 1.334939 

time= 124 ref. Vol= 149.8196 abu= 1.148048 

time= 125 ref. Vol= 149.8196 abu= 1.101325 

time= 126 ref. Vol= 43.0516 abu= 2.438934 

time= 127 ref. Vol= 43.0516 abu= 2.346022 

time= 128 ref. Vol= 86.10321 abu= 2.229882 

time= 129 ref. Vol= 86.10321 abu= 2.171812 

time= 130 ref. Vol= 30.00112 abu= 3.233213 

time= 131 ref. Vol= 10.38101 abu= 9.632974 

time= 132 ref. Vol= 20.76202 abu= 9.054996 

time= 133 ref. Vol= 20.76202 abu= 8.814172 

time= 134 ref. Vol= 20.76202 abu= 8.332523 

time= 135 ref. Vol= 5.832028 abu= 17.6611 

time= 136 ref. Vol= 11.66406 abu= 16.97523 

time= 137 ref. Vol= 11.66406 abu= 16.46083 

time= 138 ref. Vol= 11.66406 abu= 15.60349 

time= 139 ref. Vol= 4.21085 abu= 22.56076 

time= 140 ref. Vol= 8.421701 abu= 21.37335 

time= 141 ref. Vol= 8.421701 abu= 20.66091 

time= 142 ref. Vol= 8.421701 abu= 19.71098 

time= 143 ref. Vol= 2.089752 abu= 47.37404 

time= 144 ref. Vol= 1.029435 abu= 90.34086 

time= 145 ref. Vol= 2.058869 abu= 87.42664 

time= 146 ref. Vol= 2.058869 abu= 82.5696 

time= 147 ref. Vol= 0.8990695 abu= 119.0119 

time= 148 ref. Vol= 0.8990695 abu= 189.0844 

time= 149 ref. Vol= 0.8990695 abu= 180.1863 

time= 150 ref. Vol= 0.8990695 abu= 173.5127 

time= 151 ref. Vol= 0.8990695 abu= 161.2779 

time= 152 ref. Vol= 0.8990695 abu= 170.176 

time= 153 ref. Vol= 0.8990695 abu= 114.5629 

time= 154 ref. Vol= 1.798139 abu= 87.86863 

time= 155 ref. Vol= 1.798139 abu= 68.40406 

time= 156 ref. Vol= 1.798139 abu= 58.39371 

time= 157 ref. Vol= 5.394417 abu= 33.36783 

time= 158 ref. Vol= 5.394417 abu= 19.64995 

time= 159 ref. Vol= 16.18325 abu= 11.86412 

time= 160 ref. Vol= 16.18325 abu= 7.476866 

time= 161 ref. Vol= 48.54975 abu= 4.078291 

time= 162 ref. Vol= 48.54975 abu= 2.451094 

time= 163 ref. Vol= 97.0995 abu= 1.503612 

time= 164 ref. Vol= 194.199 abu= 0.8341958 

time= 165 ref. Vol= 194.199 abu= 0.5200851 

time= 166 ref. Vol= 582.597 abu= 0.288364 

time= 167 ref. Vol= 1165.194 abu= 0.1287339 

time= 168 ref. Vol= 2330.388 abu= 0.07981503 

time= 169 ref. Vol= 2330.388 abu= 0.06865809 

time= 170 ref. Vol= 2330.388 abu= 0.0622214 

time= 171 ref. Vol= 2330.388 abu= 0.05964672 

time= 172 ref. Vol= 2330.388 abu= 0.05664293 

time= 173 ref. Vol= 2330.388 abu= 0.05492648 

time= 174 ref. Vol= 2330.388 abu= 0.05321002 

time= 175 ref. Vol= 2330.388 abu= 0.04977712 

time= 176 ref. Vol= 2330.388 abu= 0.0523518 

time= 177 ref. Vol= 2330.388 abu= 0.04891889 

time= 178 ref. Vol= 2330.388 abu= 0.04763155 

time= 179 ref. Vol= 2330.388 abu= 0.04720244 

time= 180 ref. Vol= 2330.388 abu= 0.06822898 

time= 181 ref. Vol= 1049.724 abu= 0.09335784 

time= 182 ref. Vol= 294.0404 abu= 0.350292 

time= 183 ref. Vol= 588.0809 abu= 0.3298866 

time= 184 ref. Vol= 268.5301 abu= 0.3649498 

time= 185 ref. Vol= 537.0601 abu= 0.3426059 

time= 186 ref. Vol= 144.7601 abu= 0.6009941 

time= 187 ref. Vol= 289.5203 abu= 0.5802702 

time= 188 ref. Vol= 289.5203 abu= 0.5630003 

time= 189 ref. Vol= 289.5203 abu= 0.5422763 

time= 190 ref. Vol= 121.6472 abu= 0.8878134 

time= 191 ref. Vol= 121.6472 abu= 0.8384904 

time= 192 ref. Vol= 243.2944 abu= 0.780947 

time= 193 ref. Vol= 243.2944 abu= 0.7357343 

time= 194 ref. Vol= 86.89084 abu= 1.242939 

time= 195 ref. Vol= 36.6628 abu= 2.863938 

time= 196 ref. Vol= 73.32561 abu= 2.673009 

time= 197 ref. Vol= 73.32561 abu= 2.522993 

time= 198 ref. Vol= 26.28158 abu= 3.5386 

time= 199 ref. Vol= 26.28158 abu= 5.441073 

time= 200 ref. Vol= 26.28158 abu= 5.288875 

# chunk #2
demo_daphnia <- function() {
## load individual-level model (DEB and phyto, differential equations)
submodel <- daphnia_deb_phyto()

## test one outer step alone
m <- sim(submodel)
out(m)

### load individual-based model (IBM, discrete, population dynamics)
deb_ibm <- daphnia_deb_ibm()

initfunc(deb_ibm) <- function(obj) {
 p <- as.list(parms(obj))
      parms(obj)["WAM"] <- length2weight(p$SAM, p$l2w)  # weight at maturity
      parms(obj)["WON"] <- length2weight(p$SON, p$l2w)  # weight of neonates
      age <- sample(1:30, 100, repl=TRUE)
      len    <- p$SON + 2 * age/30 ## simple linear approximation
      weight <- length2weight(len, p$l2w)

      initlist <- list(X = c(x1=0, x2=0, x3=0.25),
                       Z = data.frame(
                         age        = age,
                         weight     = weight,
                         egg        = sample(1:5, 100, repl=TRUE),
                         eggage     = runif(100, min=0, max=5),
                         length     = len,
                         eggstorage = 0
                        ),
                        refVol = 10 # 10L
      )
      init(obj) <- initlist #list2vec(initlist)
      obj
}

observer(deb_ibm) <- function (statelist) {
    X <- statelist$X
    Z <- statelist$Z
    refVol <- statelist$refVol
    hist(Z$length, col="red", breaks=seq(0.5,3,0.25), ylim=c(0,150))
    c(X["x1"], X["x2"], X["x3"], abundance = nrow(Z)/refVol,
        age = mean(Z$age), weight = mean(Z$weight), eggs = mean(Z$egg),
        length = mean(Z$length), refVol = refVol)
}

#deb_ibm <- initialize(deb_ibm)

times(deb_ibm)["to"] <- 200
parms(deb_ibm)$samplesize <- c(100, 200)
parms(deb_ibm)$life.span <- 60

deb_ibm <- sim(deb_ibm)

o <- out(deb_ibm)

### plot results
par(mfrow=c(3,1))
plot(o$time, o$x1, type="l",
  main="Phytoplankton", col="red", ylim=c(0,0.5),xlab="Day",ylab="mg C / L")
lines(o$time, o$x2, col="blue")
lines(o$time, o$x3, col="green")
plot(o$time, o$abundance, type="l", main="Abundance",xlab="Day", ylab="Ind / L")
bm.mgc <- o$weight * o$abund / 1000
plot(o$time, bm.mgc , type="l", main="Biomass",xlab="Day", ylab="mg C / L")
}
demo_daphnia()

time= 1 ref. Vol= 10 abu= 10 

time= 2 ref. Vol= 20 abu= 9.8 

time= 3 ref. Vol= 20 abu= 9.3 

time= 4 ref. Vol= 20 abu= 8.8 

time= 5 ref. Vol= 20 abu= 8.1 

time= 6 ref. Vol= 2.053388 abu= 49.187 

time= 7 ref. Vol= 4.106776 abu= 46.752 

time= 8 ref. Vol= 1.740159 abu= 59.76464 

time= 9 ref. Vol= 3.480319 abu= 56.31668 

time= 10 ref. Vol= 0.6481041 abu= 160.4681 

time= 11 ref. Vol= 0.6481041 abu= 154.2962 

time= 12 ref. Vol= 0.6481041 abu= 171.2688 

time= 13 ref. Vol= 0.6481041 abu= 160.4681 

time= 14 ref. Vol= 0.6481041 abu= 197.4991 

time= 15 ref. Vol= 0.6481041 abu= 185.1555 

time= 16 ref. Vol= 1.296208 abu= 138.8666 

time= 17 ref. Vol= 1.296208 abu= 101.8355 

time= 18 ref. Vol= 2.592416 abu= 74.83366 

time= 19 ref. Vol= 2.592416 abu= 46.6746 

time= 20 ref. Vol= 5.184832 abu= 30.4735 

time= 21 ref. Vol= 5.184832 abu= 21.02286 

time= 22 ref. Vol= 15.5545 abu= 12.53657 

time= 23 ref. Vol= 15.5545 abu= 7.39336 

time= 24 ref. Vol= 31.10899 abu= 4.436016 

time= 25 ref. Vol= 62.21799 abu= 2.507313 

time= 26 ref. Vol= 124.436 abu= 1.52689 

time= 27 ref. Vol= 124.436 abu= 0.892025 

time= 28 ref. Vol= 373.3079 abu= 0.5303932 

time= 29 ref. Vol= 373.3079 abu= 0.3402017 

time= 30 ref. Vol= 746.6159 abu= 0.206264 

time= 31 ref. Vol= 746.6159 abu= 0.1366164 

time= 32 ref. Vol= 1493.232 abu= 0.1178652 

time= 33 ref. Vol= 1493.232 abu= 0.1111683 

time= 34 ref. Vol= 1493.232 abu= 0.1051411 

time= 35 ref. Vol= 1493.232 abu= 0.0984442 

time= 36 ref. Vol= 1493.232 abu= 0.09040794 

time= 37 ref. Vol= 1493.232 abu= 0.08438074 

time= 38 ref. Vol= 1493.232 abu= 0.08371105 

time= 39 ref. Vol= 1493.232 abu= 0.07567479 

time= 40 ref. Vol= 1493.232 abu= 0.06830822 

time= 41 ref. Vol= 2986.464 abu= 0.0636204 

time= 42 ref. Vol= 2986.464 abu= 0.06161133 

time= 43 ref. Vol= 669.6107 abu= 0.1403801 

time= 44 ref. Vol= 1339.221 abu= 0.1344065 

time= 45 ref. Vol= 195.5068 abu= 0.5114913 

time= 46 ref. Vol= 391.0135 abu= 0.4859167 

time= 47 ref. Vol= 186.1969 abu= 0.5639191 

time= 48 ref. Vol= 186.1969 abu= 0.5370658 

time= 49 ref. Vol= 80.25729 abu= 1.183693 

time= 50 ref. Vol= 160.5146 abu= 1.146313 

time= 51 ref. Vol= 77.9197 abu= 1.360375 

time= 52 ref. Vol= 77.9197 abu= 1.334707 

time= 53 ref. Vol= 77.9197 abu= 1.732553 

time= 54 ref. Vol= 77.9197 abu= 1.668384 

time= 55 ref. Vol= 77.9197 abu= 1.540047 

time= 56 ref. Vol= 77.9197 abu= 2.027729 

time= 57 ref. Vol= 77.9197 abu= 2.41274 

time= 58 ref. Vol= 27.92821 abu= 3.75964 

time= 59 ref. Vol= 27.92821 abu= 3.652221 

time= 60 ref. Vol= 27.92821 abu= 4.798016 

time= 61 ref. Vol= 27.92821 abu= 4.511568 

time= 62 ref. Vol= 10.01011 abu= 9.9899 

time= 63 ref. Vol= 20.02022 abu= 9.790102 

time= 64 ref. Vol= 8.483144 abu= 12.14172 

time= 65 ref. Vol= 8.483144 abu= 11.78808 

time= 66 ref. Vol= 8.483144 abu= 21.80795 

time= 67 ref. Vol= 8.483144 abu= 20.62914 

time= 68 ref. Vol= 8.483144 abu= 23.57616 

time= 69 ref. Vol= 3.70443 abu= 28.0745 

time= 70 ref. Vol= 3.70443 abu= 39.41227 

time= 71 ref. Vol= 3.70443 abu= 41.3019 

time= 72 ref. Vol= 3.70443 abu= 40.22211 

time= 73 ref. Vol= 3.70443 abu= 48.86042 

time= 74 ref. Vol= 1.722991 abu= 56.87785 

time= 75 ref. Vol= 1.722991 abu= 71.96789 

time= 76 ref. Vol= 1.722991 abu= 67.90519 

time= 77 ref. Vol= 1.722991 abu= 88.2187 

time= 78 ref. Vol= 1.722991 abu= 84.73639 

time= 79 ref. Vol= 0.640517 abu= 131.1441 

time= 80 ref. Vol= 1.281034 abu= 129.5828 

time= 81 ref. Vol= 1.281034 abu= 150.6596 

time= 82 ref. Vol= 1.281034 abu= 124.8991 

time= 83 ref. Vol= 1.281034 abu= 107.7255 

time= 84 ref. Vol= 2.562068 abu= 67.13327 

time= 85 ref. Vol= 2.562068 abu= 46.44686 

time= 86 ref. Vol= 5.124136 abu= 30.44416 

time= 87 ref. Vol= 5.124136 abu= 21.07673 

time= 88 ref. Vol= 15.37241 abu= 10.14805 

time= 89 ref. Vol= 30.74481 abu= 5.854646 

time= 90 ref. Vol= 30.74481 abu= 3.480262 

time= 91 ref. Vol= 61.48963 abu= 2.276807 

time= 92 ref. Vol= 122.9793 abu= 1.39861 

time= 93 ref. Vol= 122.9793 abu= 0.8375396 

time= 94 ref. Vol= 368.9378 abu= 0.4960186 

time= 95 ref. Vol= 368.9378 abu= 0.3903097 

time= 96 ref. Vol= 368.9378 abu= 0.3659154 

time= 97 ref. Vol= 368.9378 abu= 0.3577839 

time= 98 ref. Vol= 368.9378 abu= 0.3442315 

time= 99 ref. Vol= 368.9378 abu= 0.3333896 

time= 100 ref. Vol= 368.9378 abu= 0.3062847 

time= 101 ref. Vol= 368.9378 abu= 0.2981533 

time= 102 ref. Vol= 368.9378 abu= 0.2981533 

time= 103 ref. Vol= 368.9378 abu= 0.2873113 

time= 104 ref. Vol= 368.9378 abu= 0.2710484 

time= 105 ref. Vol= 368.9378 abu= 0.3875992 

time= 106 ref. Vol= 368.9378 abu= 0.4445194 

time= 107 ref. Vol= 368.9378 abu= 0.4499404 

time= 108 ref. Vol= 368.9378 abu= 0.425546 

time= 109 ref. Vol= 41.54705 abu= 2.238426 

time= 110 ref. Vol= 83.09409 abu= 2.118081 

time= 111 ref. Vol= 40.5337 abu= 2.368399 

time= 112 ref. Vol= 81.06741 abu= 2.09702 

time= 113 ref. Vol= 25.57331 abu= 4.692391 

time= 114 ref. Vol= 25.57331 abu= 4.457772 

time= 115 ref. Vol= 25.57331 abu= 5.083424 

time= 116 ref. Vol= 25.57331 abu= 4.770598 

time= 117 ref. Vol= 25.57331 abu= 7.625136 

time= 118 ref. Vol= 25.57331 abu= 7.507826 

time= 119 ref. Vol= 11.16739 abu= 9.491926 

time= 120 ref. Vol= 11.16739 abu= 9.13374 

time= 121 ref. Vol= 22.33477 abu= 8.417368 

time= 122 ref. Vol= 4.370797 abu= 24.70945 

time= 123 ref. Vol= 4.370797 abu= 24.70945 

time= 124 ref. Vol= 4.370797 abu= 25.85341 

time= 125 ref. Vol= 4.370797 abu= 24.93824 

time= 126 ref. Vol= 1.783999 abu= 56.61439 

time= 127 ref. Vol= 3.567998 abu= 53.8117 

time= 128 ref. Vol= 1.585777 abu= 61.79937 

time= 129 ref. Vol= 3.171553 abu= 61.16876 

time= 130 ref. Vol= 1.026393 abu= 110.0943 

time= 131 ref. Vol= 1.026393 abu= 100.3515 

time= 132 ref. Vol= 2.052785 abu= 96.45431 

time= 133 ref. Vol= 2.052785 abu= 91.09574 

time= 134 ref. Vol= 0.986916 abu= 95.2462 

time= 135 ref. Vol= 0.986916 abu= 132.7367 

time= 136 ref. Vol= 0.986916 abu= 112.4716 

time= 137 ref. Vol= 1.973832 abu= 85.11362 

time= 138 ref. Vol= 1.973832 abu= 70.42139 

time= 139 ref. Vol= 1.973832 abu= 54.20927 

time= 140 ref. Vol= 5.921496 abu= 27.86458 

time= 141 ref. Vol= 5.921496 abu= 17.0565 

time= 142 ref. Vol= 11.84299 abu= 11.31471 

time= 143 ref. Vol= 23.68598 abu= 6.332859 

time= 144 ref. Vol= 47.37197 abu= 4.137468 

time= 145 ref. Vol= 47.37197 abu= 2.70202 

time= 146 ref. Vol= 94.74394 abu= 1.519886 

time= 147 ref. Vol= 189.4879 abu= 0.9499288 

time= 148 ref. Vol= 189.4879 abu= 0.8813229 

time= 149 ref. Vol= 189.4879 abu= 0.8602133 

time= 150 ref. Vol= 189.4879 abu= 0.8391038 

time= 151 ref. Vol= 189.4879 abu= 0.7599431 

time= 152 ref. Vol= 189.4879 abu= 0.717724 

time= 153 ref. Vol= 189.4879 abu= 0.6649502 

time= 154 ref. Vol= 189.4879 abu= 0.6385633 

time= 155 ref. Vol= 189.4879 abu= 0.6121764 

time= 156 ref. Vol= 189.4879 abu= 0.5857895 

time= 157 ref. Vol= 189.4879 abu= 0.5277382 

time= 158 ref. Vol= 189.4879 abu= 0.8496586 

time= 159 ref. Vol= 189.4879 abu= 0.8179943 

time= 160 ref. Vol= 189.4879 abu= 0.7704978 

time= 161 ref. Vol= 21.03084 abu= 3.899036 

time= 162 ref. Vol= 42.06168 abu= 3.61374 

time= 163 ref. Vol= 42.06168 abu= 3.471093 

time= 164 ref. Vol= 42.06168 abu= 3.233347 

time= 165 ref. Vol= 13.48131 abu= 7.269324 

time= 166 ref. Vol= 26.96262 abu= 6.750087 

time= 167 ref. Vol= 26.96262 abu= 6.342114 

time= 168 ref. Vol= 26.96262 abu= 6.156672 

time= 169 ref. Vol= 10.29108 abu= 9.717158 

time= 170 ref. Vol= 20.58215 abu= 9.134128 

time= 171 ref. Vol= 20.58215 abu= 9.5714 

time= 172 ref. Vol= 20.58215 abu= 9.085543 

time= 173 ref. Vol= 7.73765 abu= 15.63782 

time= 174 ref. Vol= 2.803497 abu= 36.02644 

time= 175 ref. Vol= 5.606993 abu= 34.59965 

time= 176 ref. Vol= 5.606993 abu= 32.99451 

time= 177 ref. Vol= 5.606993 abu= 31.03268 

time= 178 ref. Vol= 1.704253 abu= 62.7841 

time= 179 ref. Vol= 3.408506 abu= 57.50319 

time= 180 ref. Vol= 3.408506 abu= 54.86274 

time= 181 ref. Vol= 3.408506 abu= 53.10244 

time= 182 ref. Vol= 0.8876319 abu= 112.6593 

time= 183 ref. Vol= 1.775264 abu= 108.1529 

time= 184 ref. Vol= 0.6699109 abu= 155.2445 

time= 185 ref. Vol= 1.339822 abu= 146.2881 

time= 186 ref. Vol= 0.4481009 abu= 229.8589 

time= 187 ref. Vol= 0.4481009 abu= 263.3335 

time= 188 ref. Vol= 0.8962018 abu= 218.7007 

time= 189 ref. Vol= 0.8962018 abu= 181.8787 

time= 190 ref. Vol= 0.8962018 abu= 152.8674 

time= 191 ref. Vol= 1.792404 abu= 101.5396 

time= 192 ref. Vol= 1.792404 abu= 56.34892 

time= 193 ref. Vol= 5.377211 abu= 30.12714 

time= 194 ref. Vol= 5.377211 abu= 19.15491 

time= 195 ref. Vol= 16.13163 abu= 10.78626 

time= 196 ref. Vol= 16.13163 abu= 7.438801 

time= 197 ref. Vol= 32.26327 abu= 4.21532 

time= 198 ref. Vol= 64.52653 abu= 2.820545 

time= 199 ref. Vol= 64.52653 abu= 1.54975 

time= 200 ref. Vol= 193.5796 abu= 0.9453476